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Abstract. Experiments have revealed multiple quantum oscillation frequencies 
in underdoped high temperature superconductor YBa2Cu306+5, corresponding to 
approximatdy 10% doping, which contain CuO bilaycrs in the unit cell. These unit 
cells arc further coupled along the c-axis by a tunneling matrix element. A model of the 
energy dispersion that has its roots in the previously determined electronic structure, 
combined with two-fold commensurate density waves, reveals multiple electron and 
hole pockets. To the extent that quasiparticlcs of the reconstructed Fermi surface 
have finite residues, however small, the formation of Landau levels are the cause of 
these oscillations and the bilayer splitting and warping of the electronic dispersion 
along the direction perpendicular to the CuO-planes are firm consequences. We 
explore this possibility in detail and find overall consistency with experiments. An 
important conclusion is that bilayer splitting is considerably renormalized from the 
value obtained from band structure calculations. It would be extremely interesting to 
perform these experiments for higher values of doping. We roughly expect the splitting 
of the frequencies to increase with doping, but the full picture may be more complex 
because the density wave order parameter is also expected to decrease with doping, 
vanishing around the middle of the superconducting dome. 
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1. Introduction 

The surprising quantum oscillations (QO) in both hole [DElElSlElElEllSlISlIinilll] 

and electron doped cuprates [12] have raised an important question concerning the 
ground state of high temperature superconductors [13J. Applied magnetic field between 
35 — 85T has been argued to quench the superconducting fiuctuations, at least to a 
large degree, revealing the normal state. This is not surprising in electron doped 
Nd2-xCexCu04 (NCCO), where the upper critical field Hc2 is less than 10 T, but is 
somewhat surprising in hole doped cuprates where Hc2 is extrapolated to be of order 
100 T or greater [13]. One of the striking recent findings is the observation of multiple 
quantum oscillation frequencies [3 US]. To understand QO in hole doped YBa2Cu306+5 
(YBCO) and stoichiometric YBa2Cu408 (Y124), we shall follow a reasoning based 
on broken translational symmetry with perhaps an unconventional order parameter, 
(ia;2_y2-density wave (DDW) [16]. The observed multiple frequencies should not only 
impose constraints on the theoretical models but also interpretation of experiments, as 
discussed recently [15] . Superficially similar results can be obtained within a mean field 
approximation using a spin density wave (SDW) theory, but we favour singlet DDW 
for numerous reasons discussed elsewhere [T7]. More importantly, the quasip articles 
of a singlet DDW have charge-e, spin- 1/2, and a ^^-f actor renormalized by residual 
Fermi liquid corrections. In the simplest treatment given here, we set g = 2. This 
characterization of the quasiparticles is consistent with a very recent measurement and 
its precise analysis [15] and perhaps eliminates any triplet order parameter, such as 
SDW or triplet DDW. 

Not only do the experiments involving multiple quantum oscillation frequencies 
indicate formation of Landau levels signifying finite quasiparticle residues even in 
underdoped cuprates, but also indicate coherent electron motion along the direction 
perpendicular to the CuO-plane. A bilayer Hamiltonian corresponding to YBCO was 
first written down in a paper in which an interlayer tunneling theory of superconductivity 
was proposed [TS]. This Hamiltonian was subsequently derived from a downfolding 
process in a band structure calculation [I9]. As long as the fermionic quasiparticles 
exist as excitations of the normal ground state, it is impossible to deny the existence 
of bilayer splitting, which results from the superposition of the electronic states of the 
layers within a bilayer block. For each value of momentum, there is a bonding and an 
antibonding state that are split in energy. In the original context [[18] it was argued 
that only in a superconducting state such a coherent linear superposition is possible. 
However, it is clear that the only requirement is the existence of a finite quasiparticle 
residue. An important effect discussed earlier [17] is that the phase of the DDW order 
parameter of the two layers within a bilayer block make a large difference. Even though 
the bilayer splitting can be substantial, the splitting of the Fermi surface areas for the 
out-of-phase case can be very small as compared to the in-phase case. We shall focus 
on these two alternatives amongst other considerations. 

For many years it has been argued that the normal state of high temperature 
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superconductors is incoherent, especially in the underdoped regime. Here we shall focus 
on very low temperatures, where a sharp statement can be made. The view that the 
normal state is a non-Fermi liquid appears to be at variance with the striking QO 
experiments mentioned above. We look for consistency with recent experiments [71 [T5] 
involving multiple frequencies, emphasizing of course the general aspects of a mean field 
theory. A further motivation is a measurement in a tilted magnetic field [20], where 
inconsistency of a scenario in which observed multiple frequencies arise from bilayer- 
split pockets is pointed out. The idea of probing QO with a tilted field is important, 
but our theoretical analyses are not in agreement with those presented in Ref. [20] . 

We emphasize a commensurate density wave order as the cause of Fermi surface 
reconstruction as revealed in quantum oscillation measurements, although some evidence 
for incommensuration does exist [1]. The pressing questions can hopefully be addressed 
in a simpler setting: Why should the Fermi liquid picture be valid for the normal state? 
Is the motion along the direction perpendicular to the CuO-planes (c-axis) coherent? 
Why do other experimental probes of the electronic structure paint a very different 
picture of the fermiology? In reality, no direct evidence for any kind of long range density 
wave order exists in the regime of interest to the QO measurements. Fluctuating order 
does not solve this dilemma, especially because the QO measurements require very 
large correlation lengths and nearly static order. The simplest possible explanation 
of the main aspects of the measurements call for long range order. Moreover, there 
are strong arguments from detailed fits to the measurements that the relatively high 
magnetic field is not the root of these observations [H], beyond the obvious effect of 
suppressing superconductivity. Indeed, previous NMR measurements in YBa2Cu408 
up to at least 23.2 T have shown no signatures of field induced order. Yet the QO 
measurements for this stoichiometric material are clear and unambiguous. Of course, 
NMR measurements [21] in higher fields of the order of 45T would be interesting. 
Given these larger issues and many others, it is not particularly attractive to focus 
on details such as incommensurate versus commensurate order. In any case, it was 
shown previously [T7] that within mean field theory it is quite simple to incorporate 
incommensurate order with very little change of the big picture; to go beyond mean 
field theory is quite difficult and is not particularly fruitful without a sufficiently strong 
motivation. An important point with regard to DDW is that it is hidden from most 
common probes and its existence perhaps could have gone unnoticed. 

The present manuscript is organized as follows: in Section 2 we set up the effective 
Hamiltonian and discuss bilayer splitting in Section 3. In Section 4 we discuss our results 
in a perpendicular magnetic field and in Section 5 those in a tilted field. In Section 6 
we discuss how variation of parameters provide contrasting evidence of the out-of-phase 
versus in-phase DDW order. In Sec. 7 we discuss in detail the temperature dependences 
and the oscillation magnitudes of both the magnetization and the specific heat within the 
Lifshitz-Kosevich-Luttinger formula but with Dingle factors reflecting vortex scattering 
rate in the mixed state. Section 8 contains remarks regarding unresolved puzzles. 
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We consider a tight-binding Hamiltonian, Ho, which captures correctly the bilayer 
splitting and the matrix elements between the unit cells; see Figure |2} 

2 

^0 = 5^ $^6(k)4.(k)c„,(k) -J2t±i^)ch{k)c,,{k) + h.C. 

j,k n=l j,k 

-tcJ24,j+ii^M^) + h-c. (1) 
The fermion annihilation operator c„j(k) depends on the bilayer index n within the 
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Figure 1. Bilayer structure of YBa^CusOg+i- Each unit cell contains a bilayer CuO 
block. Note that tunneling matrix elements within a bilayer unit and between the 
nearest neighbour planes between the unit cells are kept; other matrix elements are 
exponentially smaller. 

unit cell and the index j refers to the unit cell. The spin indices are suppressed and 
k = {kj., ky) is a two-dimensional vector. Only the hopping matrix element, tc-, between 
the nearest neighbor planes of two adjacent unit cells are kept, as the tunnelling matrix 
elements to further neighbour planes are considerably smaller. For simplicity t^ is 
assumed to be momentum independent, as very little is known about its precise form. 
This assumption will have little effect on our analysis. The bilyaer matrix element 
corresponding to YBa2Cu306+5 is [IH] 

i±(k) = ^ [cos(A;^a) - cos{kya)f , (2) 

where a is the in-plane lattice constant, ignoring slight orthorhombicity. Hq can be 
further simplified by the canonical transformation [22], 

c„(k, k,) = ^Y^ c„j(k)e''=^[^'^+("-i)'^]e^'<^('^^)/2^ (3) 
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which diagonahzes it in the momentum space. Note the additional phase factors, e~^'^^^''^ 
for n = 1 and e^'"'^^''"'^ for n = 2. The choice of the phase (j){kz) = —k^d preserves the 
fermion anticommutaion rules and results in an energy spectrum that is periodic in 
27r/c, which preserves the periodicity of the conventional unit cell. The 2x2 bilayer 
block is still not diagonal and must be diagonalized further to obtain the quasiparticle 
dispersion. Note that tc is a matrix element between the nearest neighbour planes of 
the two bilayer blocks and will be chosen to be an adjustable parameter. The canonical 
transformation leads to 

Ho = 5^ jj^ekcUA;, k)c„(A;„k) - [t^(k) + t,e-*'=^^]cI(A;„ k)c2(A;„ k) + h.c.\ (4) 

We use a common band structure [IS], [23] , 

= — 2t(cos k;j.a + cos kya) + 4t' cos k^a cos kya — 2t"(cos 2fc^a + cos 2kya). (5) 

with t' = 0.32t and t" = ^t', but t is chosen to be t ~ O.leV. This renormalized value of 
t, as compared to the band structure value of 0.38eV^, seems to be phenomenologically 
more appropriate in the underdoped regime of interest to us. However, the specific 
results pertaining to the ground state at T = are independent of the magnitude of t; 
even if we had chosen t = 0.38e^, the results would have been the same provided the 
remaining parameters are chosen proportionately. This is no longer be true when we 
consider the T 7^ properties discussed in Sec. 6. We shall first choose t± = 0.05t and 
tc = O.OlSt; these parameters are expected to be highly renormalized in the underdoped 
regime. Even when bilayer splitting is clearly observed in angle resolved photoemission 
spectroscopy (ARPES) in heavily overdoped Bi2Sr2CaCu208+5 (Bi2212) [24j, the actual 
magnitude of t± is severely overestimated by the band structure calculations, 300 meV, 
as opposed to the observed 88 meV. In the underdoped regime, the band structure value 
is likely to be more unreliable because of strong correlation effects. In Sec. 6 we shall 
see how the variation of t± affects the principal conclusions. 
The 2x2 Hamiltonian has the eigenvalues 

A±(fc„ k) = ek ± ^/t^c + t±i^y + 2tct±(k) cosfc.c . (6) 

It is interesting to note that with our choice of the phase of the fermion operators 
the distance between the layers in a bilayer block, d, does not appear explicitly in 
the spectrum, only implicitly in the magnitude of the hopping matrix elements. The 
above result is very different from the conventional warping of layered materials that 
contain only one electronically active plane per unit cell, which leads to a dispersion 
ek — 2tc cos kzC. It is clear that the two bilayer split bands are warped differently as a 
function of kz, as seen in Figure |2} As long as tc is nonzero, the splitting at the nodal 
locations (7r/2a, 7r/2a) is non-zero. We now fold the Brillouin zone to the reduced 
Brillouin zone (RBZ) bounded by ky ± k^ = ±7r/a, corresponding to the two- fold 
commensurate singlet DDW order parameter and augment Hq by H': 

H'= [^W^kCi(A;„k)tci(A;„k + Q) + /i.c.]+(1^2), (7) 




Figure 2. A greatly exaggerated illustration of warping of the bilayer split bands 
plotted in the extended zone — 27r < fc^c < 27r, — tt < kxO. < n, and — tt < kya < tt. 



where Q = (7r/a,7r/a) and the DDW gap is real and is given by 
Wo 

Wk = {cos k^a — cos kya). (8) 

Note that the DDW order parameters are chosen to be in-phase for the layers. To 
reproduce the experimental frequencies we require a somewhat large value of Wo = 0.85t 
within our mean field approximation; in Sec. 6 we shall consider a small variation of 
this parameter. 

The singlet DDW condensate is defined by [23] 

k^)Cn,a(M, ^z)) = « W^k 5<7',a5„',„5k',k+Q, (9) 

Note that it involves 6a',a for spin indices. This is the reason why the spin indices can 
be conveniently suppressed. This is a particle-hole condensate that breaks the following 
symmetries: translation by a lattice spacing, time reversal, parity, and a rotation by 7r/2, 
while the product of any two are preserved. The order parameter corresponds to angular 
momentum i = 2. Since there is no exchange symmetry between a particle and a hole, 
the orbital wave function does not determine the spin wave function. Therefore, there is 
also a corresponding triplet DDW, which consists of a staggered pattern of circulating 
spin currents [25], as opposed to a staggered pattern circulating charge currents. In 
the present work, we shall consider only the singlet DDW order and make only brief 
remarks regarding the triplet DDW at the very end. The staggering is determined by 
the wave vector Q. 
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3. Bilayer splitting 



The combined Hamiltonian Hq + H' can be written in terms of tlie four component 
spinor \l/"l'(k, /c^) = {c\{k,kz),c\{k + Cl,kz),cl{k,kz),cl(k + Q,kz)}, suppressing once 
again the spin indices, which is irrelevant for a singlet DDW order parameter. In terms 
of this spinor the combined Hamiltonian is 

M= ^Hk,A;,)A^(k,A;^) (10) 

where 

/ ek iWk -t^ik) - t^e-'^''^ \ 

-t^ik) - tce'""^^ Ck iWk 

\ -t^{]^) - t.e'ck. -iW^ ek+Q y 

Note that the DDW order parameters in the two 2x2 diagonal blocks are in phase. The 
in-phase DDW order parameter corresponds to "ferromagnetically" aligned staggered 
circulating currents in the layers within a bilayer block. The four eigenvalues of the 
matrix A are 



.fll^ 



and 
A2±(k) 



± 



'-^-^^^ + Wl - ^tl + t+(k)2 + 2tM^) COS K 



,(12) 



Ck + Ck+Q 




Ck — Ck+Q 



+ ^tl + t^ikY + 2t,t^ik) COS kzc (13) 



For a particle-hole condensate, when measured from the chemical potential, the spectra 
are 

i?,^±(k) = At±(k)-/i, (14) 
i?2\(k) = A^(k)-/i. (15) 

because both e(k) and e(k + Q) are equally shifted by /i. If, on the other hand, 
the staggered circulating currents are "antiferromagnetically" aligned within a bilayer 
block [TTj, that is, iW^ is replaced by —iW]^ in the lower 2x2 diagonal block [out-of- 
phase), the corresponding eigenvalues are 



A?±(k) 

and 

A^±(k) 



Ck + Ck+Q 



Ck + Ck+Q 



Ck — Ck+Q 



Ck — ek+Q 



^tl + t+(k)2 + 2tct±{k) cos A;^c 



+ t+(k)2 + 2tct±(k) cos fc^c 



1/2 



(16) 



1/2 



.(17) 
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Figure 3. Bilayer splitting of the Fermi surfaces for the in-phase DDW order 
parameter. The cut is at fc^ = where the spUtting is maximal. For clarity the 
contours are plotted in the extended zone. 

Once again, measured from fi , we have 

^i±(k) = At±(k)-/i, (18) 

EL(k) = A^,i(k)-/i. (19) 

The contour plots for the Fermi surfaces corresponding to Xl^ and A|_|_ for = are 
shown in Figure [3j It is clear that while the electron pockets are observably split, the 
splitting of the hole pockets is much smaller. The chemical potential, n = — 0.78t, was 
adjusted to yield approximately 10.3% hole doping. For identical set of parameters, the 
splitting for the out-of-phase eigenvalues, A^_|_ and Af-t, is considerably smaller, as shown 
in Figure |4j It would be incorrect, however, to infer that the splitting is exactly zero; 
see Table [2] below. Note that the absolute value of t does not change the frequencies 
because ep Oi',at,at',at",at±,atc,aWo) = aep Oi]t,t' ,t" ,t±,tc,Wo), as long as we also 
let fi —J- a/i. 

4. Magnetic field perpendicular to the CuO-plane: Onsager relation 

Quantum oscillation frequencies can be obtained from the extremal areas, Alep), of the 
Fermi surface perpendicular to the applied magnetic field . The Onsager relation for 
the frequency F is 

F=^^M. (20) 

Of course, this formula presupposes that the quasiclassical approximation is valid and 
there are no significant magnetic breakdown effects. 
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Figure 4. Bilayer splitting of the Fermi surfaces for the out-of-phase DDW order 
parameter. The cut is at kz = 0. For clarity the contours are plotted in the extended 
zone. 



Given the electronic structure, the doping dependence can be obtained from noting 
that there are two hole pockets within the RBZ and one electron pocket. These are 
further split by the bilayer coupling and warped by the dependence. Taking into 
account two spin directions, the doping fraction of a given electron pocket corresponding 
to the bilayer bonding band 6, x\ is 

:-— / dK / dky / dk, 9 {fi - e\k, k,)) . (21) 

i^-^r Jo Jo J-n/c 

There is an identical expression for the antibonding contribution Xg. Similarly the two 
hole pockets contribute an amount x^ given by 



4 



< = ^77^, / dk^ / dky / dk, 9 (e^(k, K) - fi) , (22) 

Jo Jo J-tt/c 



(2vr; 



with an identical antibonding contribution x^. The total hole doping per CuO-plane is 
then 

Xh = 2'-"^^ + ~ ~ ^e)- (23) 

The frequencies for the in-phase order parameter are given in Table [T] The 
parameters were chosen, but not particularly optimized, to be similar to the observed 
frequencies [7] 540 ± 15T, 630 ± 40T, 450 ± 15T and 1130 ± 20T. Out of 4 theoretically 
predicted frequencies corresponding to the electron pocket only 3 are observed. The 
fourth observed frequency at 1130T could correspond to the hole pocket that is split 
very little. Alternately, it may also be a harmonic. It has been puzzle for some 
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Table 1. Bilayer split frequencies for the in-phase DDW order. Here t± — 0.05i, 
tc — O.OlSt and doping is approximately ^ 10.3%. The band parameters are given 
in the text. The electron pocket is labelled as e-pocket and the hole pocket as the 
h-pocket. 

e-pocket {k^ = 0) e-pocket {k^ = tt/c) h-pocket {kz = 0) h-pocket (fc^ = tt/c) 

711 T 659 T 1051 T 1032 T 

480 T 534 T 997 T 1015 T 



Table 2. Bilayer split frequencies for the out-of~phase DDW order. The parameters 
are the same as in Table [TJ 



e-pocket {k^ = 0) 


e-pocket (fc^ = tt/c) 


h-pocket {kz = 0) 


h-pocket {kz = n/c) 


617 T 


609 T 


1044 T 


1031 T 


573 T 


585 T 


1002 T 


1016 T 



time [271 EHl EHl HH |30l [3T] as to why the hole pocket frequencies have such weak 
or non-existent signatures in quantum oscillation measurements. 

In contrast, the out-of-phase frequencies (Table [2]) do not resemble the experimental 
observations [7]. We provide an alternate picture in Sec. 6 based on the experiment 
in Ref. [IS]. Within the mean field approximation adopted here, it is not possible to 
distinguish between the in-phase and the out-of-phase cases [I7| as far as the electronic 
energy is concerned. For this one would need a detailed microscopic Hamiltonian. This 
is outside the scope of the present investigation. We therefore rely on experiments to 
distinguish between the two cases. 



5. Tilted magnetic field 

In this section we calculate the effect of tilted magnetic field on quantum oscillations [32] . 
In Figure [s] we show a cut of the Fermi surface with the plane ky = —n/a. The 
intersection with the plane A is given by 

ep {k^,ky,k^tainp) = fi (24) 

If the Fermi surface does not depend kz, the area Ao will be constant for all planes 
perpendicular to k^ = 0, and the area in the plane A will be given by Aa = Ao/ cos 9? 
with a constant value of Aq- However, given the dependence on kz, the area is 

/ — 7r/fi p—n/a 
dkx / dky6 [ep {kx, ky, k^ tarn/)) — yu] , (25) 
-■n/a J — ixja 

which can be computed numerically. The above result corresponds to /c^ = 0. More 
generally, when the plane O is situated at an arbitrary value of kz, we get 

dkx / dkyO [ep {k^, ky, kz + k^ tan (p) — fi] . (26) 

■Tv/a J—w/a 
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Figure 5. A cut of the Fermi surface of a poclcet by a plane fcy = — vr/a. B is tiie 
applied magnetic tilted at an angle tp. 



Note that Equation 26 is valid for angles < (fr. 

vr/c a 
vr/a c 



such that 



tan Lfr. 



(27) 



Beyond this maximum angle, there are discontinuous jumps, and we do not attempt to 
treat this case. For hole pockets the frequencies are summarized in Figure |6} Similarly 
for the electron pockets the frequencies are shown in Figure [7} Note that the vertical 
scales are different in Figures |6] and [7j It is also interesting to note how the frequencies 
are sequentially split as we first turn on t± and then tc, which is shown in Figure |8] for 
the magnetic field in the direction kz- 

In Figure [9] we illustrate the warping along for an electron pocket. We represent 
a cut for kx = that yields kya around 2.64. The vertical line at the center corresponds 
to a model where t± = tc = 0. The two vertical lines at the left and the right of the figure 
show the splitting when t± ^ 0. The line on the left corresponds to the outer pocket 
and the line on the right to the inner pocket. As we turn on tc, the warping is seen as 
two curved lines. It is clear that the warping has opposite sense for the outer and the 
inner pockets. The bilayer splitting can be seen from the displacement of the left line 
by 0.0334 from the central line, while the right line is displaced by 0.0379 in the oposite 
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Figure 6. Hole pocket frequencies as a function of the tilt angle (p. 
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Figure 7. Electron pocket frequencies as a function of the tilt angle ip. 



direction. The splitting induced by t± is therefore not symmetric. To calculate the 
warping we can compute distance between the lines at = 0, obtaining 0.00937 for the 
outer pocket and 0.0119 for the inner one. At kz = ±7r/c, 0.00958 is the displacement 
for the outer pocket and 0.01158 is the displacement for the inner one. These numbers 
encode two important facts: first, the warping is different for the inner and the outer 
pockets, and, second, it cannot be modeled with a simple cosine dependence. 



6. Variation of parameters 

Here we vary the parameters to see how the results change. The focus is the difference 
between the out-of-phase and the in-phase DDW order parameters. We have already 
seen that there is a qualitative distinction between them. However, given the recent 
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Figure 8. The hierarchy of frequency splitting for in-phase DDW order parameter 
with the magnetic field normal to the CuO-plane, as we sequentially turn on t± and 
tc, not to scale. 



measurements [15], we would like to see if one or the other can be made more consistent 
with these experiments. We stress that the phenomenological nature of our work 
precludes us from fitting parameters with certainty, nor is it our intention. We only 
look for some qualitative insights. However, since in this section we shall be computing 
the oscillatory part of the thermodynamic potential, as a function of temperature and 
magnetic field, not just the frequencies, a good estimate of the leading tight-binding 
matrix element t is necessary for materials relevant for quantum oscillation experiments. 
Since there are no reliable ARPES for YBCO, the next best we can do is to rely on 
the recent tight-binding fit to the measured ARPES in Y124 [33], a system in which 
good quantum oscillations have been observed. Except for t, the ratios of the remaining 
band parameters to t are not very different from the band structure results given below 
Equation |4j Thus, we simply take over the value of t determined from ARPES in Y124, 
which is t = 0.154 eV (the average of the fit to the bonding and the antibonding bands). 
Additionally, we would like to see if one can tolerate a much larger value of bilayer 
matrix element, t±, as compared to the earlier section and still find consistency with 
experiments. We shall see that this is indeed possible, but only for the out-of-phase 
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Figure 9. Warping of the bilayer split electron pockets. 



Table 3. Bilayer split frequencies for the out-of-phase DDW order parameter. Here 
t± = 0.12t and doping is « 10.7%. The band structure parameters, t' and t" , are 
unchanged, t' = 0.32t and t" = 0.5i, but Wq = OM, where t = 0.154 eV. 



e-pocket {kz — 0) 


e-pocket (fc^ — n/c) 


h-pocket {kz = 0) 


h-pocket {kz = n/c) 


538 T 


535 T 


1034 T 


1015 T 


461 T 


474 T 


975 T 


993 T 



DDW order parameter. 

In this section we keep all the band structure parameters fixed, including tc, but 
more than double the bilayer matrix element to t_i_ = 0.12t = 0.0185 eV, resulting in 
a splitting of 37 meV, which is reasonable compared to the overdoped Bi2212, where 
it is measured to be 88 meV; one expects renormalization with underdoping. To keep 
the doping level more or less fixed (~ 10.7%), we set /x = — 0.775t and Wq = 0.9t. The 
resulting oscillation frequencies are shown in Table |3} and the Fermi surfaces at A;^ = 
are plotted in Fig. [TT] The two groups of electron pocket frequencies are close to each 
other and so are the two groups of hole pocket frequencies despite much larger bilayer 
splitting. The warping of the outer electron pocket is only 3 T and that of the inner 
pocket is 13 T. It is even possible to tolerate larger t±, but we have not explored it 



further. It is again useful to examine the frequency diagram. This is shown in Fig. [TO 
and is quite different from Fig. [8] 

In contrast, for the same set of parameters, as above, the in-phase DDW results in 
frequencies that are no longer close to the recent experiments [T3] as shown in Table |4| 
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Figure 10. The hierarchy of frequency splitting for out-of-phase DDW order 
parameter with the magnetic field normal to the CuO-plane, as we sequentially turn 
on t± and tc, not to scale. 



Table 4. Bilayer split frequencies for the in-phase DDW order parameter. Here 
t± ~ 0.1 2t and doping is « 10.7%. The band structure parameters are unchanged, 



t' = O.m and t" = 


0.5i, but Wo 


= OM, where t = 0.154 eV. 


Compare with Table [sj 


e-pocket {k^ = 0) 


e-pocket (fc^ 


= tt/c) h-pocket (fc^ = 0) 


h-pocket {kz — tt/c) 


750 T 


699 T 


1046 T 


1022 T 


250 T 


307 T 


965 T 


988 T 



The warping of the outer electron pocket is 51 T and that of the inner electron pocket 
is 57 T. 

6.1. Cyclotron masses and the second derivatives of the extremal areas 

In the following section we shall need the cyclotron masses and the second derivatives 
of the extremal areas. These are calculated numerically and are sumarized in Table |5] 
and Table [6j The second derivatives of the extremal areas with respect to kz are more 
difficult to calculate. We fit the areas near kz = and /c^. = vr by a fourth order 
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Figure 11. Bilayer splitting for the out-of-phase DDW order parameter. The cut is 
at kz = 0. For clarity the contours are plotted in the extended zone. Here t± — 0.12t 
and doping is ~ 10.7%. The band structure parameters, t' and t", are unchanged, 
t' = 0.32t and t" = OM, but Wq = 0.9t, where t = 0.154 eV. 

Table 5. Cyclotron masses in units of the free electron mass for the out-of-phase 
DDW order parameter. 



e-pocket {k^ = 0) 


e-pockct {kz = tt/c) 


h-pocket {kz = 0) 


h-pockct (kz = tt/c) 


1.67 


1.64 


0.99 


0.97 


1.47 


1.49 


0.93 


0.95 



polynomial with only even terms (odd terms are zero within numerical precision) 

A^Ao + -^ihcr + ^ihcr (28) 

A^Ao + ^ iihc) + ^ liKc) - vr]^ (29) 

It is interesting to note that while the cyclotron masses depend on the in-plane hopping 
matrix element t, the second derivatives of the extremal areas are independent of t. 

7. Oscillation amplitudes of specific heat and magnetization 

Within Fermi liquid theory Luttinger [M] has shown that the thermodynamic potential 
is given by (/3 = l/ksT) 
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Table 6. The fitting coefficients of tlie extremal areas. 









A2 


A4 


h-pocket (fcz 


= 0) 


1.45542 


-0.00804 


-0.00051 


h-pocket {kz 


-0) 


1.37238 


0.00776 


-0.00488 


e-pocket {kz 


-0) 


0.75696 


-0.00119 


-0.00305 


e-pocket {kz 


-0) 


0.64928 


0.00902 


-0.02421 


h-pocket {kz 


= tt/c) 


1.43000 


0.04116 


-1.16330 


h-pocket {kz 


= tt/c) 


1.39761 


-0.04102 


1.18022 


e-pocket {kz 


= tt/c) 


0.75315 


0.00269 


-0.01303 


e-pocket {kz 


= tt/c) 


0.66771 


-0.00862 


0.00637 



where {Er} constitute the spectra of elementary excitations behaving hke independent 
particles in a magnetic field, including Fermi liquid corrections; r denotes the collection 
of quantum numbers: the Landau level n, k^, and the spin a. The spectra {Er} in a 
crystalline solid in high magnetic fields are of course not easy to calculate, especially if 
we have to include bilayer sphtting and the DDW order discussed above, but a rigorous 
answer can be given within an asymptotic expansion. Luttinger has shown that the 
problem maps onto to that solved by Lifshitz and Kosevich (LK) |35] in which the 
thermodynamic potential depends on the extremal areas of closed orbits, the derivative 
of the areas with respect to energy at the chemical potential, and the second derivative 
of the extremal areas with respect to k^. The beauty of this approach is that it is 
not necessary to know Er explicitly. Thus, even given the complexity of the present 
problem, the procedure to calculate the oscillatory part of the thermodynamic potential 
is straightforward. As with all asymptotic expansions, the validity of the procedure 
far surpasses what we may naively perceive to be the regime of validity. Thus the LK 
formula has stood the test of time, especially with Luttinger's Fermi liquid corrections. 

For simplicity, in this section we shall consider magnetic field only in the c-direction. 
Taking into account only the fundamental frequencies, Fi, the oscillatory part of Q is 

-1/2 



Q 

V 



oc 



dkl 



cos 



2TiFi vr 

± - 

H 4 



cos vr 



ni 



m 



(31) 



The phase ±7r/4 correspond to positive or negative sign of the second derivative of the 
extremal area with respect to kz- The sum is over all extremal surfaces, and is the 
cyclotron effective mass given by (m is free electron mass) 



2tx 



(32) 



92 s, 



dm 



is the second derivative of the area of the Fermi surface with respect to kz. The 



argument of the function 



sinh A,- 



(33) 
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The cyclotron frequencies are given by u^^ = 
heat is then 



(34) 

eH/m*c. The oscillatory part of the specific 
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^1/2 
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(35) 



where ifj" (A) is the second derivative of ip (A)with respect to A. Similarly, the leading 
oscillatory term of the magnetization is 

-1/2 



dpi 



^(A. 



cos 









± - 


cos vr — 
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4_ 


\ m J 



.(36) 



These results need to be supplemented by the Dingle factors that damp the oscillations 
due to scattering from defects or vortices in the vortex liquid state or both. We expect 
that the total scattering rate to be given by the combination of defect and vortex 
scattering rates 

h 



_h h 

Moreover, these scattering rates must depend on the particular extremal area, 
consideration. The calculation of the Dingle factors 



(37) 
under 

(38) 



especially in the mixed phase including disorder, with coexisting fluctuating ci-wave 
superconducting order parameter and DDW, is a daunting task. Previously, we have 
shown rigorously that almost any form of conventional disorder due to defects in a pure 
DDW state suppresses the electron pockets more than the hole pockets. For the vortex 
scattering rate, however, an approximate treatment based on a paper by Stephen 
led to an interesting prediction relating the Dingle factors of electron and hole pockets 
(not including bilayer splitting) in the commensurate case, which is 



h 



V2 



m 
m* 



3/2 



h 



UJcTv 



(39) 



where m is a characteristic scale having the dimension of mass corresponding to the 
massless nodal fermions of the DDW (!), and m* is the cyclotron mass corresponding 
to the electron pocket as defined above (note that the notations are different here from 
Ref. [21]), which in turn is very close to the band mass defined by expanding around 
the bottom of the electron pocket. Although the precise numerical relation is difficult 
to control, it is reasonable to set = T^ei with a = 1.5 — 4.5 for phenomenological 
purposes; we had estimated this parameter earlier to be 4.4 
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7.1. Specific heat and magnetization 

With the frequencies given in Table [3} the cyclotron masses in Table |4| and the second 
derivatives of the areas in Table [5| we can compute the oscillatory parts of the specific 
heat and the magnetization provided we can make reasonable estimates of the Dingle 
factors. The Dingle factor for electrons is a bit more controlled because the band mass 
obtained obtained by expanding around the antinodal points is quite consistent with 
the computed cyclotron masses. Assuming that samples have negligible disorder, we 
shall estimate the scattering rate of the electrons to be given by the vortex scattering 
rate, which, following an analysis of Stephen was found to be [TTl 122] 



1-7^ A ^ (40) 



Jv) e ^ V J V \fJ'\hUc 

where Aq is the magnitude of the T = superconducting gap, which we set to be 
PS lOmeV for the relevant doping range. The cyclotron frequency Uc = eH/m*c, with 
m* given in Table |4j With the present set of parameters, and with the average value of 
m*, we find that 

1 



5 X 10'^ s"\ (41) 



where we used, as a typical case, if = 40 T and Hc2 ~ 100 T. We believe that this 
gives the correct order of magnitude; for the earlier set of parameters we estimated it to 
be 3 X 10^^ s^^ [17J. The Dingle factors of the holes are more complex p!7l [29] because 
it has to be estimated taking into account the nodal fermions for DDW, but it is not 
unreasonable to assume a ~ 2 in the relation Vh = 'D'^, uniformly for all electron and 
hole pockets. 

The computed specific heat at four representative temperatures are shown in 



Fig. 12, It is interesting to note that there is a vr-phase shift from high to low 



temperatures. The same results are visualized in a 3D-plot in Fig. 13 As a function 
of temperature and magnetic field the oscillations go through a node, which is also 
visible in Fig. 12 The reason for this is the factor ip'^X) in the formula for the specific 
heat. The Fourier transform of the oscillations in 1/H on the other hand shows a more 
complex structure for specific heat as shown in Fig. 14, which, however, is very sensitive 
to the Dingle factor, the range oi 1/H over which the Fourier transform is performed, 
and the windowing technique. The results shown here uses no windowing technique, and 
the range of the magnetic field is 1/60 < 1/H < 1/20. The non-monotonic behavior of 



the Fourier transform in Fig. 14 can be understood by glancing at Fig. [13] Because of 
the aforementioned node, the transition from 6K to 4K lowers the amplitude. At 2K 
the amplitude recovers again and then finally decreases again at IK. Note that only 
one dominant frequency is seen. 

Similarly, we also plot the oscillations of the magnetization as a function of 1/H 
in Fig. 15, but it is difficult to detect multiple frequencies with naked eyes. Even in 
the Fourier transform over a range 1/60 < 1/H < 1/20, shown in Fig 16, the multiple 
electron pocket frequencies known to be present in the formula are not resolved. The 



Multiple quantum oscillations 



20 




Figure 12. Oscillatory part of the specific heat as a function of the magnetic field 
at representative temperatures. Note the phase shift by tt as the the temperature is 
lowered. The vertical scale is in arbitrary units. 




6 



Figure 13. Oscillatory part of the specific heat as a function of temperature and 
the inverse of the magnetic field. Note the presence of nodes for intermediate values 
magnetic field and temperature. The vertical scale is in arbitrary units 



Fourier transform is now monotonic as a function of temperature unlike the results for 
the specific heat. The arrow indicates weak, unresolved hole pocket frequencies around 
1000 T; see Fig. 10 The two electron pocket frequencies at 535 T and 538 T strongly 
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Figure 14. Absolute value of the Fourier transform, /, in arbitrary units of the 
oscillations of the specific heat as a function of the frequency F in units of Tesla. No 
windowing was performed, and the field range was 1/60 < 1/H < 1/20. 




Figure 15. Oscillatory part of the magnetization as a function of temperature and 
the inverse of the magnetic field. There are no nodes at intermediate temperatures as 
in the case of specific heat. The vertical scale is in arbitrary units 

overlap and lead to a large amplitude. The spin interference factor cos(7rm*/m) play 
an important role in the respective weights of the various frequencies. 
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Figure 16. Absolute value of the Fourier transform, /, in arbitrary units of the 
oscillations of the magnetization as a function of the frequency F in units of Tesla 
with no windowing and the field range of 1/60 < 1/H < 1/20. 

8. Conclusion 

We have argued that bilayer sphtting and warping of the electronic dispersion in kz are 
necessary consequences of a reconstructed Fermi liquid normal state, and measurements 
in a tilted magnetic field can be useful in probing the frequency spectra of quantum 
oscillations. The small value of the warping is intimately connected to the large ratio of 
the in-plane to c-axis resistivities. However, the magnitude of bilayer splitting necessary 
to produce overall consistency with experiments for the in-phase DDW order is strongly 
renormalized (~ 10 meV) from the band structure value (~ 300 meV). Note that the 
distance between the layers is only 3.25A, similar to the in-plane lattice constant. In 
contrast, with the out-of-phase DDW order a larger value of bilayer splitting (~ 37 meV) 
can be tolerated. This is an important consequence of the out-of-phase DDW order. 
Although strong electronic correlations in the underdoped regime must be responsible 
for such renomalized parameters, a convincing explanation is missing despite many 
speculations, especially because the effective mass is only about twice the free electron 
mass. It would be interesting to carry out these QO measurements for larger hole- 
doping for which we generally expect the splitting to increase, unless some other effects 
involving the decrease of the magnitude of the order parameter intervenes. It is worth 
emphasizing once again that even in heavily overdoped Bi2212, the renormalization of 
the observed bilayer splitting, 88 meV, in ARPES, as compared to the band structure 
value of 300 meV, is still not understood. 

The calculations presented here can be easily extended within a mean field theory 
to SDW and incommensurate order along the lines discussed elsewhere [T7j. A more 
illuminating exercise is to compare and contrast quantum oscillations in hole and 
electron doped cuprates [37]. The likely differences in the upper critical fields lead 
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to important physical differences. Further work along this direction is in progress. The 
triplet DDW [25] at the simplest mean field level produces results similar to SDW, 
which is also a triplet order parameter, but with orbital angular momentum zero. Such 
triplet order parameters are necessary to explain the experiments |6j involving the non- 
existence of spin zeros in QO. However, more recent experiments [12] have revealed spin 
zeros and have concluded that quasiparticles behave like charge-e, spin-1/2 fermions with 
a (jf-factor consistent with 2.2. This is strongly indicative of a singlet order parameter, 
but not a triplet order in the particle-hole channel, such as SDW or triplet DDW [38] . 

Although we have obtained consistency with experiments using Fermi liquid theory, 
it is not certain that non-Fermi liquid aspects should be ignored, at least insofar as 
underdoped YBCO is concerned. Convincing explanation of the lack of the hole pocket 
frequencies required by the Luttinger sum rule [SH HOI SI] and the inconsistency with 
Fermi arcs observed in ARPES, albeit in zero magnetic field, are intriguing. We know of 
one example, the z/ = 1/2 quantum Hall effect, which despite being a non- Fermi liquid 
has a phenomenology similar to a Fermi liquid in many respects [42j. The situation in 
NCCO is clearer [37], however. We hope that our work will shed light on these exciting 
set of experimental developments. 
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